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ABSTRACT 



Context. In photo-dissociation regions (PDRs), polycyclic aromatic hydrocarbons (PAHs) may be produced by evaporation of very 
small grains (VSGs) by the impinging UV radiation field from a nearby star. 

Aims. We quantitatively investigate the transition zone between evaporating very small grains (eVSGs) and PAHs in several PDRs. 
Methods. We studied the relative contribution of PAHs and eVSGs to the mid-IR emission in a wide range of excitation conditions. 
We fitted the observed mid-IR emission of PDRs by using a set of template band emission spectra of PAHs, eVSGs, and gas lines. 
The fitting tool PAHTAT (PAH Toulouse Astronomical Templates) is made available to the community as an IDL routine. From the 
results of the fit, we derived the fraction of carbon /evso locked in eVSGs and compared it to the intensity of the local UV radiation 
field. 

Results. We show a clear decrease of /evsG with increasing intensity of the local UV radiation field, which supports the scenario of 
photo-destruction of eVSGs. Conversely, this dependence can be used to quantify the intensity of the UV radiation field for diff'erent 
PDRs, including unresolved ones. 

Conclusions. PAHTAT can be used to trace the intensity of the local UV radiation field in regions where eVSGs evaporate, which 
correspond to relatively dense (hh = [100, 10^] cm~^) and UV irradiated PDRs (Go = [100,5 x 10"^]) where H2 emits in rotational 
lines. 

Key words. ISM: Photon-dominated region (PDR); Astrochemistry; ISM: dust, extinction; ISM: lines and bands; Infrared: ISM; 
ISM: molecules 
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1. Introduction 

The mid-infrared (mid-IR) spectrum of photo-dissociation re- 
gions (PDRs) is dominated by the aromatic infrared bands 
(AIBs) attributed to the emission of polycyclic aromatic hydro- 
carbons (PAHs, Leger & Puget 1984; Allamandola et al. 1985; 
Joblin & Tielens 2011). Very small grains (VSGs) are another 
class of mid-IR emitters if they are small enough to be stochas- 
tically heated (Desert et al. 1990). The VSGs are generally con- 
sidered to be carbonaceous (Desert et al. 1986) but Li & Draine 
(2001b) concluded that the presence of very small silicate grains 
is not excluded and can involve up to 10% of interstellar Si. 
Cesarsky et al. (2000) suggested a link between PAHs and VSGs, 
in which the latter are processed by the UV radiation field or by 
shocks giving birth to free PAHs. 

Using blind- source separation methods, Rapacioli et al. 
(2005) and Berne et al. (2007) extracted three independent spec- 
tra from the mid-IR emission, two molecular PAH-type spectra 
attributed to neutral and ionised PAH populations, respectively 
and one consisting of both continuum and broad band emissions, 
which was attributed to some VSGs. The relative spatial distri- 
bution between PAHs and these VSGs additionally supported 
the idea that free flying PAH molecules are produced by photo- 
evaporation of VSGs under UV irradiation at the border of PDRs 
(Rapacioli et al. 2005; Berne et al. 2007). Although they are 
called VSGs, this population is not strictly that included in the 
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Desert et al. (1990) model: it is better described by the PAH 
grains of Li & Draine (2001a). In the following, we will refer to 
them as evaporating VSG (eVSGs). 

Using the template spectra derived by Rapacioli et al. (2005) 
and Berne et al. (2007), Joblin et al. (2008) and Berne et al. 
(2009b) were able to fit the mid-IR spectra of objects in dif- 
ferent evolutionary stages, such as planetary nebulae and proto- 
planetary disks. A spectrum attributed to large ionised PAHs 
(PAH^) was introduced to fit highly irradiated objects. In Berne 
et al. (2009b) and Berne et al. (2009a), the relative intensities of 
the different mid-IR components were related to physical quan- 
tities such as the spectral type of the central illuminating star for 
proto-planetary disks, or the ionisation parameters and radiation 
field for the Monoceros R2 compact HII region. While it is clear, 
in all the above studies, that the evolution of eVSGs into PAHs 
is related to the intensity of the UV field, this has never been 
quantified. 

In this paper, we search for a quantitative link between the 
intensity of the UV radiation field and the chemical evolution 
of the AIB carriers in terms of processing of eVSGs into PAHs. 
Section 2 describes the PDR sample. Section 3 presents the fit- 
ting tool we have developed to decompose mid-IR spectra into 
chemical populations, which takes into account dust extinction. 
We apply this method to our PDR sample and present the re- 
sults in Sect. 4. Section 5 describes the relation between the lo- 
cal physical conditions and the evolution of the AIB carriers and 
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illustrates the potential of our method for characterising the UV 
field in a variety of PDR environments. 

2. The PDR sample 

We analysed mid-IR observations of several PDRs that span a 
wide range of UV irradiation conditions. The majority of the 
observations discussed here were obtained with the short-low 
(SL) module of the Infrared Spectrograph (IRS) (Houck et al. 
2004) onboard Spitzer (Werner et al. 2004) in spectral mapping 
mode. Data reduction was performed with the CUBISM soft- 
ware (Smith et al. 2007a) and consisted of cube assembling, 
calibration, flux correction for extended sources and bad pixel 
removal. For the objects that were not observed with IRS, we 
used ISOCAM data available in the ISO data archive as highly 
processed data products (Boulanger et al. 2005). 

2.1. PDRs with spatially resolved mid-IR populations 

Our sample contains several PDRs (Table 1) where the emission 
from the PAH and eVSG populations (see Sect. 3) can be spa- 
tially isolated. These PDRs are referred to as PDRs with spatially 
resolved mid-IR populations. It is possible to study the variations 
of the mid-IR spectral properties in these sources (Sect. 3) and 
link them to the local physical properties such as the local den- 
sity and UV field (Sect. 5) as a function of the depth inside the 
PDR. To estimate the intensity of the UV radiation field at the 
PDR front, we used modelled stellar spectra from Kurucz (1993) 
and applied a dilution factor based on their projected distance to 
the PDR front. All parameters used to calculate the intensity of 
the UV radiation field at the PDR fronts are reported in Table 
1. In the following, we will express the UV field intensity Go 
in terms of the Habing field, which corresponds to 1.6 x 10"^ 
erg cm"^ s"^ when integrated between 91.2 and 240 nm (Habing 
1968). The geometrical model assumed for these PDRs is dis- 
cussed in Appendix A. 2. 

p-oph filament Assuming typical spectral properties for the B2 
star HD 147889 and using geometrical dilution to the projected 
distance of 0.4 pc between the star and the PDR leads to an in- 
tensity of the UV radiation field of Go ~ 520. 

NGC 2023 For NGC 2023 an extinction of Ay = 1.1 mag was 
applied, reflecting absorbing dust around the illuminating star 
HD 37903 (Compiegne et al. 2008). Using a distance of 0.38 
and 0.13 pc for the north and south PDRs, respectively yields 
an estimate of the UV radiation field intensity of Go ~ 400 and 
Go ~ 4000. These estimates are consistent with previous values 
obtained in the far-IR study by Burton et al. (1998). 

NGC 7023 The estimate of the UV radiation field intensity at 
the PDR front in NGC 7023 NW is more uncertain. Chokshi 
et al. (1988) estimated Go ~ 2.4 x 10^ through far-IR observa- 
tions and modelling, but values as high as Go ~ 10"^ have been 
proposed to explain the observed H2 ortho-to-para ratio (Fuente 
et al. 1999) towards this PDR. This value agrees with the value 
determined by assuming typical stellar properties for the corre- 
sponding spectral type of the illuminating binary star (Alecian 
et al. 2008) and geometrical dilution to the PDR front. However, 
extinction between the star surface and the PDR may attenuate 
the radiation field, yielding a lower value at the PDR front. In this 
work, an extinction correction of Ay = 1.5 was derived from the 



lUE spectrum measured on the star and was applied on the cor- 
responding Kurucz spectrum. Assuming that the same extinction 
factor applies between the star and the PDR leads to an estimate 
for the NW PDR of Go ~ 2.6 x 10^ which is consistent with the 
value of Chokshi et al. (1988). 

2.2. PDRs with mixed mid-IR populations 

In the previous section we have defined the PDRs with spa- 
tially resolved mid-IR populations as those where the PAHs and 
eVSGs can be spatially isolated. The PDRs of our sample that do 
not follow this criterion are referred to as PDRs with mixed mid- 
IR populations. The fact that in these sources the emission from 
PAHs and eVSGs cannot be isolated may result from diff'erent 
eff'ects, e.g. from the small spatial extent of the PDR, a lower 
signal-to-noise ratio in the data, specific geometrical properties, 
etc. For each of these sources we analysed an averaged spectrum 
towards the peak of the mid-IR emission. The sample of PDRs 
with mixed mid-IR populations comprises a wide range of irra- 
diation conditions, from Go ~ 100 for the Horsehead Nebula to 
Go ~ 4 X 10^ for the Orion Bar. For these well- studied objects, 
a reliable estimate of the intensity of the UV radiation field can 
be found in the literature. 



3. A spectral fitting tool for studying the chemical 
evolution of AIB carriers: PAHTAT 

Our motivation here is to provide a tool that allows us to decom- 
pose observed mid-IR (5.5-14yum) spectra into a limited set of 
components that can provide information on the chemical nature 
of the emitters. We note that there already exists an AIB fitting 
tool called PAHFIT, proposed by Smith et al. (2007b), which 
decomposes mid-IR spectra using a collection of lorentzians 
and black body functions with a large number of parameters. 
While PAHFIT provides high-quality fits that are useful e.g. to 
extract the integrated intensity from individual bands, our aim 
is to decompose the observations using band template spectra 
for PAHs and eVSGs, minor additional unattributed PAH fea- 
tures, gas lines, and a function to describe the continuum. The 
corresponding IDL code PAHTAT (PAH Toulouse Astronomical 
Templates) is available for download^. 

3.1. Band emission 

Applying blind- signal separation methods to ISO and Spitzer 
spectro-imagery data of several PDRs Rapacioli et al. (2005) 
and Berne et al. (2007) extracted from the observations a set 
of three elementary mid-IR spectra that can reproduce the AIB 
emission. Joblin et al. (2008) and Berne et al. (2009b) found that 
a fourth spectral component was needed to explain the observed 
spectra in highly irradiated environments. PAHTAT is based on 
these four astronomical template band spectra (Fig. 1). Here, we 
briefly present the spectral properties of this set of templates and 
their chemical assignment. Details on the spectral characteris- 
tics of the templates (e.g. relative band strength, position, width 
etc.), how they were obtained as well as the populations they rep- 
resent can be found in Joblin et al. (2008). These templates can 
be downloaded on the PAHTAT webpage. The first component 
is the neutral PAH (PAH^) spectrum, which is mainly charac- 
terised by a high C-H (11.2 yum) over C-C (7.7 yum) band ra- 
tio. The second component is the PAH cation (PAH^) spectrum. 



^ http : //userpages . irap . omp . eu/~ c joblin/PAHTAT. 



2 



R Pilleri et aL: Evaporating Very Small Grains as tracers of the UV field in PDRs 
Table 1. Observational parameters for PDRs with spatially resolved mid-IR populations. 



Object 


Star 


Distance 
LPCJ' 


Spectral 
Type 


Kurucz 
Spectrum * [K] 


Radius * 


^front 
I J 


/^front 


Instrument 


NGC 7023 NW 
NGC 7023 S 
NGC 7023 E 


HD 200775 


430 


B3Ve — B5^^'"'^^ 


2 X 15000^^) 


10 


42 
55 
155 


2600 
1500 
250 


Spitzer IRS-SL 
ISO-CAM cvf 
ISO-CAM cvf 


NGC 2023 N 
NGC 2023 S 


HD 37903 


470 


BL5V^ 


24000 


6 


164 
67 


400 
4000 


Spitzer IRS-SL 
Spitzer IRS-SL 


p-Oph filament 


HD 147889 


118 


B2III — B2V 


22000 


5 


610 


520 


ISO-CAM cvf 



* Derived from spectral type. 

(a) van den Ancker et al. (1997) — (b) Racine (1968) — (c) Finkenzeller (1985) — (d) Witt et al. (2006) — (e) Alecian et al. (2008) 
(f) Diplas & Savage (1994) 



Table 2. The Spitzer IRS observational parameters and the assumed physical parameters for PDRs with mixed mid-IR populations. 



Object 


Position 

[0^2000, <^2000] 


Aperture 


Go 


[cm-3] 


Ref. 


Horsehead 


(05:40:53.8; -02:28:00) 




100 


2x 10^ 


Habart et al. (2005) 


Ced 201 


(22:13:25; 70:15:03) 


22" X 22" 


300 


4x 102 


Young Owl et al. (2002) 


IC 63 


(00:59:01; 60:53:19) 


30" X 40" 


1100 


1 X 10^ 


Gerin et al. (2003) 


Parsamyan 18 S 


(06:59:41; -07:46:45) 


9" X 8" 


3500 


1 X 10^ 


Ryder et al. (1998) 


Parsamyan 18 N 


(06:59:41; -07:46:12) 


4" X 8" 


5000 


1 X 10^ 


Ryder et al. (1998) 


Orion Bar 


(05:35:21.4; -05:25:15) 


4" X 8" 


4x 10^ 


2x 10^ 


Tauber et al. (1994) 



which displays on the contrary a low 1 1.2 fim over 7.7 fim band 
ratio. The third PAH component is associated to large (~ 100 C 
atoms) ionised PAHs, referred to as PAH^, whereas the PAH^ 
and PAH^ populations would be dominated by smaller species. 
Its spectrum shares the same characteristics as that of PAH^ but 
has a "7.7" yum band that is shifted by 0.15 yum to the red com- 
pared to PAH^ (for a discussion of the reason for this shift, see 
Joblin et al. 2008). The fourth component is the eVSG spectrum, 
whose carriers are still discussed and have typical sizes of ~ 500 
C atoms (Rapacioli et al. 2005). PAH clusters (Rapacioli et al. 
2006) or complexes of PAHs with Fe atoms (Simon & Joblin 
2009) have been proposed. eVSGs are thought to carry contin- 
uum emission, but here we have treated band and continuum 
emission separately. The eVSG band spectrum is made of broad 
features at 7.7 and 1 1.2 yum and has no 8.6 yum band (Fig. 1). We 
improved this spectrum as compared to Joblin et al. (2008) by 
adding a Gaussian plateau centred at 12.6yum with a FWHM of 
1.76yum and a relative intensity compared to the 7.7 fim band of 
0.36. This plateau was introduced to better fit the extracted VSG 
component of Rapacioli et al. (2005) and Berne et al. (2007). 
The continuum model is described in the following section. 



3.2. The continuum 

In Rapacioli et al. (2005) and Berne et al. (2007) it was shown 
that some continuum emission is spatially associated with the 
band emission of eVSGs. In mild UV field PDRs, the continuum 
emission is observed to be linear up to ~ 20yum (Berne et al. 
2007; Compiegne et al. 2008). In high UV field PDRs, such as 
the Orion Bar, the shape of the continuum is observed to change 
due to a contribution from big grains at thermal equilibrium. For 
most of our sources we found that the 5.5-14 jim continuum can 
be fitted by a single slope. For high UV field PDRs a bilinear 
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10 12 14 

Wavelength (|am) 



Fig. 1. Template spectra used in the fitting procedure. The eVSG 
spectrum from Joblin et al. (2008) was slightly improved to bet- 
ter fit the observed spectrum by adding a plateau in the 12-14 yum 
range. Spectra were normalised to their integrated intensity be- 
tween 5.5 and 14 jim. 



continuum is used. The choice of the continuum (single slope or 
bilinear) can be set by the user in PAHTAT. 
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3.3. Minor PAH features 

Three additional features are present in the mid-IR spectra of 
most PDRs, at 6.7, 11.55 and 13.5 jim. These features are best 
represented by three lorentzians at 6.7, 11.55 and 13.5 jim with 
a FWHM of 0.2, 0.4 and 0.25 yum, respectively. These features 
also appear in the PAH^ and PAH^ components of Rapacioli 
et al. (2005). However, they are much fainter than the principal 
PAH features and were not included in the band template spectra. 
Their intensities are left as free parameters in the fit. However, 
these features are minor and their introduction does not change 
the results significantly. 

3.4. Gas lines 

Several gas lines fall in the 5.5-14yum range, such as the H2 S(2) 
to S(5) rotational lines, and fine structure lines from ionised gas 
such as [Aril], [Nell] and [SIV]. Some of these lines may be 
blended with PAH/eVSG features at the spectral resolution of 
ISOCAM-cvf and IRS-SL. Therefore, these features cannot be 
simply subtracted but they rather need to be treated simultane- 
ously with the PAH and eVSG bands. In PAHTAT, we assumed 
Gaussian profiles with a fixed central wavelength and a FWHM 
given by the average spectral resolution R of the instrument (R 
= 45 for ISOCAM, and R = 80 for IRS). 




Wavelength (|Lim) 

Fig. 2. Eff'ect of the extinction correction on a linear continuum 
and on a synthetic PDR spectrum composed of continuum, band 
emission, and gas lines. 



3.5. Extinction correction 

Finally, our model takes into account extinction along the line 
of sight by dust within the PDR, which is caused mainly by sili- 
cate grains. Assuming that the emitting and absorbing materials 
in the considered column of material are fully mixed^ (see, for 
instance Disney et al. 1989; Smith et al. 2007b) leads to a correc- 
tion factor (1 - ^"^^)/t^. The optical depth at each wavelength 
is calculated using the relation 



(1) 



where Cext(^) is the extinction cross-section per H nucleon, cal- 
culated by Weingartner & Draine (2001) for 7?v = 5.5. 

3.6. Summary of the model 

In summary, the modelled mid-IR emission /model is written as 

/model ~ (/bands ~^ /continuum ~^ /gas ~^ /mbands) 7^ 77 • (^) 
^ ^ Cext ^Vh 

The adjusted parameters are /bands ^ which represents the intensity 
of the band template spectra of PAH^'^'^ and eVSGs, /continuum, 
which is the (bi-)linear continuum emission, /gas, which con- 
sists of the entire mid-IR H2 and ionised gas lines, and /mbands 
to account for the minor band features at 6.7, 1 1.5 and 13.5 yum. 
The hydrogen column density A^h is left as a free parameter in 
the fit. Hereafter, A^h will be expressed, for simplicity, in mag- 
nitudes of visual extinction along the line of sight, defined as 
Ay = A^h/(1.8 X 10^^ cm"^) mag. We tested the eff'ect of using 
diff'erent extinction curves, corresponding to diff'erent values of 
Ry, and found a diff'erence of up to 20% in the estimated column 
density but no significant eff'ect on the population weights when 
using Ry = 3.1 instead of Ry = 5.5. The extinction correction 
can significantly modify the shape and intensity of the mid-IR 



^ PAHTAT gives the possibility to also use a foreground extinction 
correction e~'^^, although we assumed a mixed extinction model here. 



spectrum. As an example. Fig. 2 shows the eff'ect of the extinc- 
tion correction to a linear continuum and a synthetic spectrum 
for diff'erent values of Ay. This shows an example of how diff'er- 
ent (i.e. non-linear) continua can be obtained using just a linear 
continuum and an extinction correction. On the other hand, to 
obtain significant information on the column density, this must 
be high enough to modify the spectrum significantly. This hap- 
pens at an Ay > 3, which we defined as our threshold. 

3.7. lo, the emission radiated by AIB carriers 

The total emission radiated by AIB carriers in both aromatic 
bands and mid-IR continuum of eVSGs (cf. Sect. 3.2) is given 
by 



/O - /b; 



I p 

bands ^continuum ~ bands bands continuum 



eVSG 



/, 



PAH 



+ /; 



eVSG _^ ^eVSG 



(3) 



where /^^^s intensity of all the PAH populations (PAH^, 

PAH^ and PAH^). if^^^ and ^^^.^ are the contribution of 

^ bands contmuum 

eVSGs to the bands and the continuum, respectively. While 
If^l^ is a direct output of PAHTAT, /^^^.^ has to be com- 

bands ^ ' contmuum 

puted separately. Considering that in moderate UV field regions 
(Go < 4000) emission from big grains is negligible at /I < 14 yum, 
we considered that for these PDRs /^^^^^^^ = /continuum- Hence 
using the data in this region, we can calibrate a relation between 
^c^ntinuum KanS ' Figure 3 shows the pixel-to-pixel correla- 
tion of the two extracted components I^^^ and /continuum for all 
resolved PDRs (where Go < 4000). The dispersion of this cor- 
relation suggests that the band and continuum emissions are not 
strictly related and that the band-to-continuum ratio can be af- 
fected by factors that still have to be understood, such as chemi- 
cal composition and size distribution. To quantify the continuum 
intensity we therefore used the limits given by the cloud of points 
in Fig. 3, which gives /^^^.^ = (1 ± 0.6) x I?^^^. The contri- 

^ ' ^ contmuum ^ ^ bands 

bution of eVSGs to emission at longer wavelengths (A> 14 fim) 
is neglected based on the results by Draine & Li (2007), who 
showed that it is not significant for PAH grains of size < lOA 
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(~ 1200 C atoms), which includes the size of eVSGs (Rapacioli 
et al. 2005). 



1 •ID- 



S' 10" 



Jfi 6-10- 



4«10" 



2-10- 



X NGC 7023 NW 
)K NGC 7023 E 
o NGC 7023 S 
□ NGC 2023 N 

NGC 2023 S 

Oph-W 



leVSG //leVSG . 
'Continui/m' 'Bands 




4«10"^ 6-10" 

lS[Wm-sr-] 



MO" 



Fig. 3. Pixel-to-pixel correlation of the two components I^^^ 
^c^ntinuum rcsolvcd PDRs of Table 1 . For each PDR, a 



mask was applied to show only the brightest pixels in / 
the pixels towards the PDR front (see also Fig. 10). 



bands 
eVSG' 



I.e. 



After defining the total eVSG emission, we can use in the 
following the ratio of this value to the total AIB emission as a 
tracer of the evolution of eVSGs in diff'erent environments: 



/e 



eVSG 



^eVSG _^ ^eVSG 
bands continuum 



(2 ± 0.6) X 



pYSG 
bands 



(4) 



which can be derived directly from PAHTAT. 
4. Results of PAHTAT 

We used PAHTAT to fit the emission spectra of the PDRs with 
spatially resolved (Table 1) and mixed (Table 2) mid-IR popula- 
tions. 

4.1. PDRs with spatially resolved mid-IR populations 

We applied PAHTAT for each spatial pixel of the spectro- 
imagery cubes of the resolved PDRs listed in Table 1 . This yields 
the spatial distributions of the diff'erent emission components, as 
well as the column density in each pixel. Figure 4 shows the 
maps obtained for NGC 7023 NW. Maps for the other PDRs are 
displayed in Fig. 9. The spatial distributions of the PAH^, PAH^ 
and eVSG components are found to be close to those reported 
in previous studies (Rapacioli et al. 2005; Berne et al. 2007; 
Rosenberg et al. 2011). The comparison between the distribu- 
tions of Ay and the AIB emission is consistent with the interpre- 
tation that the mid-IR emission is dominated by the excited edge 
of a dense structure (cloud or filament). The column densities 
determined with this method are in line with those derived from 
other tracers (see for example Habart et al. 2003; Gerin et al. 



1998) and agree in absolute values within a factor of 4. The spa- 
tial patterns of the Ay maps seem indeed to be associated with 
variations in the physical or geometrical properties. For instance, 
the peak of column density observed in NGC 7023 NW perfectly 
matches the high-density region traced by CS emission (Pilleri 
et al., in prep.). 

4.2. PDRs with mixed mid-IR populations 

Figure 5 shows the observed spectra of the PDRs with mixed 
mid-IR populations and their fits, whose main results are re- 
ported in Table 3. The quality of the fit can be evaluated us- 
ing the signal-to-noise ratio (S/R), defined as the ratio between 
the power of the fit divided by the power of the residuals (see 
Table 3). All S/Rs are above 100 (20 dB), which attests to the 
good quality of the fit. The extinction we derive is compati- 
ble with previous estimates that used molecular tracers or (sub- 
)millimeter dust emission (see, for example, Tauber et al. 1994; 
Kemper et al. 1999; Habart et al. 2005; Karr et al. 2005). The 
weights of the diff'erent populations vary in the PDRs, reflect- 
ing a difl'erence in chemical abundances. PAH^ were suggested 
to be present only in highly excited regions (Joblin et al. 2008; 
Berne et al. 2009a), and not in mild UV-field reflection nebulae 
(RapacioH et al. 2005; Berne et al. 2007). Indeed, PAH^ repre- 
sent a minor fraction of the AIB emission (<10%) except for the 
nucleus of the starburst galaxy M82. 



5. Linking the chemistry with the physical 
conditions 

Rapacioli et al. (2005) and Berne et al. (2007) have suggested 
that eVSGs might be PAH clusters that evaporate into free-flying 
PAH molecules when exposed to UV radiation in PDRs. In the 
following, we assume that all carbon contained in eVSGs is 
transferred into PAHs, i.e. that the number of C atoms per unit 
volume contained in AIB carriers, hq = n^^^ -\- n^^, is con- 
served. We have shown in Joblin et al. (2010) that there are two 
independent ways to derive the abundance of C in the mid-IR 
carriers, either by analysing the intensity of the mid-IR emission 
or by using the total column density derived from the fit. We 
have checked that the ratio of these two quantities varies by less 
than 10% in the cut of NGC7023 NW shown in Fig. 10, exclud- 
ing strong variation in the abundance of carbon at the transition 
eVSG/PAH. 

In Appendix A.l we show that the emissivity per carbon 
atom of eVSGs and PAHs can be considered as similar. This 
means that /evsG directly represents the fraction of carbon atoms 
locked up in eVSGs compared to all carbon atoms in the AIB 
carriers. Therefore, we can quantify the evolution of the abun- 
dance of eVSGs as a function of the physical conditions (Go and 
the gas density ^h) by using the quantity /evsG, which can be 
directly derived from the observations (cf. Equation 4). 



5.1. The relation between Go and /evsG 

For each of the spatially resolved PDRs, we focused on one 
cut that extends perpendicularly from the star's position to the 
PDR and a few arcsecs past the PDR front (see Figs. 4 and 
9). Qualitatively, the variation of /evsG, calculated using Eq. 4, 
along each of the cuts shows that the eVSG abundance increases 
from the PDR front into the inner layers and reaches its max- 
imum just behind the front. The carbon content in eVSGs de- 
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21^0l"'36'0 30'0 
RA (2000) 



Fig. 4. Maps of the different components obtained using R\HTAT to fit the observed mid-IR spectra of the NGC 7023 NW IRS 
spectral cube. Integrated intensities in the 5.5-14 yum range of the total emission (left map) and of the different band template 
components (central panels) are displayed. The map of the column density expressed in Ay is also shown, as well as the flux in the 
H2 S(3) line (right maps). The black arrow shows the direction of the spatial cut studied in this work (Fig. 10). Contours of the total 
emission map are given in steps of 3 x 10"^ W m"^ sr"^ and are reported for reference in all the maps. The color intensity scales are 
in units of lO'^W m ^ sr \ while the column density is expressed in magnitude of visual extinction. Results for the other resolved 
PDRs can be found in Fig. 9. 




Wavelength (|am) Wavelength (|am) 

Fig. 5. Mid-IR spectra (black) obtained with IRS and their fit with PAHTAT (red) for the PDRs listed in Table 2. All quantities 
derived from the fit are reported in Table 3. The dashed line represents the continuum determined by the fit, corrected for the 
extinction along the line of sight. 



creases significantly in the more exposed part of the PDRs, con- 
firming the photo-evaporation scenario. 



To quantify this process, we need an estimate of the value 
of Go at various depths of the PDR along this cut, which is 
not known a priori. This was performed using the approach 
presented in Habart et al. (2003) for the Oph-W filament. This 
method relies on the modelling of the AIB emission using a few 



physical and geometrical assumptions^. This yields the values 
of Go and nu at various depths in the PDR (upper panels in Fig. 
10), which can be compared with the variation of /evsG (lowest 
panels in Fig. 10). 

In Fig. 6, we report the values of /evsG as a function of the in- 
tensity of the UV radiation field Go for all regions. For each cut, 
the number of points that are reported depends on the spatial res- 



^ For the sake of completeness and to facilitate reading, the descrip- 
tion of this method is reported in Appendix A. 
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Table 3. Results of PAHTAT for the PDR with mixed mid-IR populations. 



Source 


Integrated Intensities [W m ^ sr 






PAH and eVSG fractions [%] 




Results 






/obs 


/bands 


/continuum 


/gas 




PAH^ 


PAH+ 


PAH^ 


eVSG 


S/R(^) 


f (c) 
JeVSG 


Av 


Horsehead 


8.3 X 10-6 


4.1 X 10-6 


3.7 X 10-6 


3.9 X 10- 


-7 


32 


2 


1 


64 


238 


0.78 


6.5 


Ced 201 


2.4 X 10-5 


1.6 X 10-5 


7.6 X 10-6 


1.2 X 10" 


-7 


21 


21 


11 


45 


281 


0.63 


< 3 


IC 63 


7.2 X 10-6 


4.2 X 10-6 


2.8 X 10-6 


2.1 X 10- 


-7 


54 








45 


111 


0.63 


< 3 


P18S 


4.8 X 10-5 


3.2 X 10-5 


1.5 X 10-6 


3.5 X 10- 


-7 


48 


18 


10 


22 


303 


0.37 


10.0 


P18N 


5.1 X 10-5 


3.4 X 10-5 


1.5 X 10-5 


5.0 X 10- 


-7 


57 


13 


11 


17 


173 


0.30 


21.2 


Orion Bar 


9.9 X 10-3 


3.9 X 10-3 


5.8 X 10-3 


4.8 X 10- 


-5 


59 


32 





8 


161 


0.15 


8.9 



Notes. Relative to /bands- ^/^ = / ^lodeM^I I (^modeiW - lobsWf Fraction of C atoms locked in eVSGs. See Eq. 4 for details. 
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Fig. 6. Fraction of carbon atoms locked in eVSGs relative to the total carbon in the AIB carriers (derived with PAHTAT) as a 
function of the intensity of the local UV radiation field. The error bars result from the propagation of the uncertainty on the intensity 
of the eVSG continuum (see Sect. 3.7). 



olution of the observations. In each spatially resolved PDR, we 
excluded the low-density (^h < 5 x 10^ cm"^) layers that lie in 
the most exposed region of the PDR (up to an Ay ~ 0. 1 - 0.3). In 
these regions, eVSGs are underabundant because the photoevap- 
oration process is expected to be much faster than reformation. 
Hence their contribution to the mid-IR emission is minor and the 
derivation of their abundance is subject to large uncertainties. 
Instead, we concentrated on the PAH-eVSG transition, which 
arises at Ay ~ 0.5 and is spatially coincident with the H2 pure 
rotational emission. This corresponds to regions in which fresh 
matter from dense regions is exposed to UV photons. Finally, 
we also added one point for each of the PDR with mixed mid-IR 
populations, for which we relied on the values of Go and nu that 
are commonly used in the literature (Table 2). 

Figure 6 shows a clear decrease of the fraction of carbon in 
eVSGs with increasing intensity of the UV radiation field. The 
trend can be fitted by the expression 

/evsG = (-0.23 ± 0.02)logio(Go) + (1.21 ± 0.05). (5) 

This implies that /evsG can be used to trace the (local) in- 
tensity of the UV radiation field. At low UV field intensities 
(Go < 100), the mid-IR fits are less reliable because the AIB 



emission is fainter and sometimes barely detected, therefore the 
/eVSG vs Go relationship cannot be established. Similarly, the 
extrapolation of the relation between /evso and Go to more 
intense UV fields (e.g. above 5 x 10^, which is our limit in 
this paper) is also difficult because in these environments, the 
photo-processing of the eVSG and PAH populations is extreme 
and leads to the difficulty to observe the more fragile eVSG 
species, and likely breaks the assumption on the conservation 
of nc in the AIB carriers. We conclude that the procedure that 
is described in this paper to derive Go is reliable in the range 
100 < Go < 5 X 10"^. The gas densities that have been explored 
span the range 10^ - 10^ cm"^. 

5.2. Correlation with other tracers 

Figure 7 displays the correlation between /evsG and the local 
density for the PDRs with both mixed and spatially resolved 
mid-IR populations, showing that these two quantities are un- 
correlated. Similarly, no correlation appears between /evso and 
the ratio Go/nu, despite the link of this ratio with the fraction of 
ionised PAHs by Fleming et al. (2010) in their detailed study of 
NGC 7023 NW, NGC 2023 N, and NGC 2023 S. This proves 
that the evaporation of the eVSG population is indeed directly 
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Fig. 7. Relation between the local density and the /evsG value for 
the PDRs with spatially resolved ( Table 1) and mixed (Table 2) 
mid-IR populations. Error bars are calculated as in Fig. 6. 



linked to the intensity of the UV radiation field through a mech- 
anism that diff'ers significantly from those driving the ionizasion 
balance. This strongly suggests that the reformation of eVSGs is 
not efficient in these environments. Rapacioli et al. (2006) have 
indeed shown that the photo-evaporation of PAH clusters in the 
NGC 7023 NW PDR is much faster than their reformation by 
collisions with PAHs. 

5.3. Use of PAHTAT to derive Go or other parameters 

In principle, PAHTAT can be used to fit any mid-IR spectrum 
that shows AIB emission to provide /evsG and therefore the lo- 
cal value of Gq. However, when using this tool, one must bear 
in mind for which frame this relationship is applicable. PAHTAT 
was calibrated in relatively dense and highly UV-irradiated en- 
vironments that are found in star-forming regions. This implies 
1) that the exciting and evaporating radiation field results from 
massive stars so that the radiation field is dominated by UV pho- 
tons, and 2) the presence of dense regions where eVSGs are 
abundant and start to evaporate where the UV irradiation pen- 
etrates. In other types of regions, the value of Go derived from 
PAHTAT may not be reUable. Still, PAHTAT can be used to iden- 
tify unusual spectral properties of AIBs as related to peculiar 
properties of the source. Vega et al. (2010) for instance used the 
PAHTAT templates and identified an exceptionally intense 11.3 
Jim features in early-type galaxies. They have interpreted this as 
a population of large neutral PAHs, which are not included in 
our templates. 

The mid-IR emission of star-forming galaxies mainly arises 
from unresolved dense and highly UV-irradiated PDRs in which 
PAHTAT could be used to derive an average radiation field. To 
test the reliability of the value of Gq derived with PAHTAT 
on unresolved PDRs, we used the fitting tool to determine a 
value of /evsG for the average spectra calculated over more ex- 
tended fields of view (up to 1 arcmin) in our template object, 
NGC 7023 NW. This can be compared with < /evso the mean 
value obtained by averaging results on the individual pixels con- 
tained in the same region. The two values agree well for regions 
close to the eVSG peak. For regions close to the PDR front, the 
value of < /evsG > is higher but by less than 0.05. This yields an 
overestimate of the mean UV radiation field by less than 50%. 
Although this error is not directly applicable to other PDRs, it 
shows that the method can be used to derive an eff'ective UV ra- 
diation field intensity in unresolved sources within a factor ~ 2. 
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Fig. 8. Mid-IR spectra (black line) of the nucleus and the outflow 
of M82 and their fits (red line) obtained with PAHTAT. The re- 
sults of the fit suggest that the UV radiation field in the nucleus 
is slightly higher, as well as the column density of the gas. 

As an example, we applied our fitting procedure to the IRS 
spectra of the nucleus and outflow of the starburst galaxy M82. 
The fit results are shown in Fig. 8. The M82 IRS spectral map 
has been studied previously by Beirao et al. (2008). They found 
that the magnitudes of visual extinction were ranging between 
Ay = to 50 and were likely responsible for most of the spectral 
variations in the AIB features. The Ay values derived from our 
decomposition on each pixel are consistent with these values. 
The results of the fit applied to the mean spectra of the nucleus 
and the outflow of M82 are still consistent with the results from 
Beirao et al. (2008), because we derive an estimate of the extinc- 
tion along the line of sight of Ay ~ 21 for the nucleus and ~ 1 1 
for the outflow. The values of /eysG derived from the mid-IR fit 
correspond to a UV radiation field of Go = 3^^ x 10^ for the 
nucleus and Go = (2 ± 1) x 10^ for the outflow. Previous studies 
using chemical modelling of the nucleus of M82 by Fuente et al. 
(2008) have predicted a value of Go ~ 1 x 10"^, whereas those 
using far-IR observations (Colbert et al. 1999; Kaufman et al. 
1999) report lower values (10^-^ < Go < 10^-^). These agree 
well with our estimates. 



6. Conclusion 

We have presented a new method for analysing the mid-IR spec- 
tra of PDRs using a basis of PAH and eVSG template spectra. 
We have shown that with a few, physically meaningful parame- 
ters it is possible to obtain a very good fit of the mid-IR spectra 
of PDRs that span a wide range of physical conditions. The fit re- 
turns the relative contribution of the PAH and eVSG populations 
to the mid-IR emission, the intensity of the gas lines (especially 
the rotational lines of H2), and the magnitudes of visual dust ex- 
tinction along the line of sight. 

We derived the abundance of C atoms locked in eVSGs rel- 
ative to the total abundance of carbon in the AIB carriers. This 
value strongly correlates with the local UV field intensity and 
can therefore be used to quantify this important physical param- 
eter. Our study has provided new insights into the properties of 
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eVSGs, supporting a chemical scenario in which these species 
are destroyed by UV photons to produce free-flying R\Hs. PAH 
clusters and PAH-Fe complexes have been proposed as plausi- 
ble carriers for these eVSGs and detailed studies on the survival 
of these species in PDRs have been performed (Rapacioli et al. 
2006; Simon & Joblin 2009). One important conclusion of our 
observational study is that eVSGs are destroyed in PDRs faster 
than they can be reformed. This suggests that they are formed 
in the denser parts of molecular clouds, although this scenario 
will be difficult to probe since in these regions mid-IR emission 
cannot be triggered and therefore eVSG fingerprints cannot be 
observed. Another intriguing aspect is the continuum emission 
that seems to be associated to eVSGs. We showed that the in- 
tensity of this continuum is quite variable relative to the band 
intensity, likely reflecting a chemical or size distribution effect. 
Solving this question opens new perspectives for studies on lab- 
oratory analogues of eVSGs. 

Finally, our fitting code PAHTAT (PAH Toulouse 
Astronomical Templates) is made publicly available to the 
community as an IDL package, together with the template 
spectra used in this work. This procedure combined with the 
high spatial resolution and sensitivity of forthcoming missions 
such as JWST and SPICA will allow us to determine the spatial 
variation of the UV radiation field in more distant galaxies and 
protoplanetary disks. 

Appendix A: Mid-IR emission model 

In this section, we describe how we derived the local physical 
parameters, i.e. the gas density (^h) and the UV radiation field 
(Go) from the local mid-IR emission of the AIB carriers. 



Energy balance 

The UV absorption cross-section of PAHs was measured to be 
proportional to the number of involved carbons (Joblin et al. 
1992). This has also been shown on a larger sample of species 
using theoretical calculations'^ (Malloci et al. 2004, 2007). Using 
a typical value for cr^^, one can estimate that screening atoms 
inside a grain by surface atoms is negligible for grains with lesser 
than -5x10"^ carbon atoms. As mentioned in Sect. 5, eVSGs 
considered in these work are relatively small, containing up to 
about 1200 carbon atoms (~ 10 A), and we therefore consid- 
ered that the absorption cross sections per C atom of eVSGs and 



PAHs are similar: cr^J^^ 



^PAH 



crc- 



The power absorbed by the AIB carriers can be expressed as 



abs 



J91 



nccrc(A)F(A)dA, 



(A.1) 



2nm 



where F(A) is the ffux of the radiation field from the nearby star. 
Since ionised PAHs and large neutral PAHs (and hence eVSGs) 
can absorb significantly in the visible (Salama et al. 1996), we 
assume Aj^ax = 1000 nm. 

As discussed by different authors (Boulanger & Perault 
1988; Habart et al. 2003), the power emitted by the AIB carriers 
can be expressed as a function of Gq: 



6 X X Go 



(A.2) 



where Go denotes the UV radiation field in Habing units. 
Considering that all species absorb with the same ctq and ap- 
plying the conservation of energy, an emissivity per C atom of 



Star 




Observer 



Fig. A.l. Spherical shell geometry used to model the PDRs with 
spatially resolved mid-IR populations (Table 1). The r and Ay 
axis refers to the local position inside the PDR, along the star- 
PDR axis, counted from the star and from the PDR front, respec- 
tively. The d axis refers to the position of the line of sight on the 
sky, counted from the star position. 



the AIB carriers can be defined as e^^^ = 
sivity 6 can be conveniently expressed as 



^PAH 



€c. The emis- 



where hq is the number of C atoms per cubic centimeter that is 
locked up in the AIB carriers. Writing the energy balance Pem = 
Pabs leads to 



-j^ MOOOnm 
J 91. Inm 



nc/nn 



crc(A)F(A)dA. 



(A.3) 



http://astrochemistry.ca.astro.it/database/ 



We computed e for a variety of PAHs of different sizes 
and charges using the values of crc available from the on-line 
database of Malloci et al. (2007). The results scatter within 10%. 
Assuming a value of nc/nu = 4.2 x 10"^ (Draine 2003, for 
Ry = 5.5) we obtain 6 ~ 5 x 10"^^ WH"\ consistent with pre- 
vious estimates by Habart et al. (2001) and Habart et al. (2003). 
Eq. A.3 shows that the value of 6 depends a priori on the spec- 
tral shape (hardness) of the radiation field. To estimate the vari- 
ation of 6 with the hardness of the radiation field, we calculated 
it with the appropriate Kurucz spectra for the stars of Table 1 . 
Our calculations show that for the objects studied in this paper, 
6 marginally depends on the hardness of the radiation field. 

A.2. Determining of the structure of tlie PDRs witli spatially 
resolved mid-IR populations 

The concepts presented in the previous section can be used to 
study the physical conditions in PDRs with spatially resolved 
mid-IR populations (Table 1). However, to directly relate the ob- 
servations to the local UV radiation field and gas density, some 
knowledge on the geometry of the PDR is needed. We assumed 
a spherical shell model as in Habart et al. (2003, see also figure 
A.l), in which the cloud is divided into successive layers of in- 
creasing density. The star is placed at the centre of the cavity, at 
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Table A.l. Parameters for the gas density profiles in the PDRs with spatially resolved mid-IR populations, extracted from our 
emission model as explained in Appendix A.2. References: (a) Gerin et al. (1998) (b) Fuente et al. (1995) (c) Habart et al. (2003) 



Object 


<^front 


plateau 


< 


a 


do 


L 




[pc] 


[cm-3] 


[10^ cm-3] 




[pc] 


[pc] 


NGC 7023 NW 


0.087 


2 X lO^^^) 


1.1 


2.0 


0.016 


0.13 


NGC 7023 S 


0.11 


2 X 10^^^^ 


3.4 


2.0 


0.054 


0.10 


NGC 7023 E 


0.32 


3 X lO^^^) 


1.4 


2.7 


0.091 


0.27 


NGC 2023 N 


0.38 


2 X 104(^) 


2.4 


2.0 


0.069 


0.27 


NGC 2023 S 


0.13 


2 X 104(^) 


3.6 


2.0 


0.033 


0.07 


Oph-W 


0.40 


4 X 10^^^^ 


2.6 


2.5 


0.051 


0.28 



a distance Jfront that is determined by the projected distance on 
the sky (cf. Table 1). The total intensity radiated by AIB carriers 
as a function of the position on the sky is obtained by integrating 
Eq. A.2 along a given line of sight, i.e.: 



1 

~ 47r Jo 



6 X ^h(^) X Go(r) J/, 



where z is the projected distance from the front, r is the distance 
from the star, and L is the thickness of the PDR perpendicu- 
larly to the plane of the sky. The local radiation field inside the 
PDR depends on the impinging radiation field at the PDR surface 
and on the gas density profile that determines its attenuation. It 
evolves as 



Go(r) = G 



front 



X e 



(A.4) 



where (Tuv is the dust extinction cross-section at 1000 A(1.5 x 
10"^^ cm^ H"\ Weingartner & Draine 2001). The density pro- 
file is modelled between the front position at r = Jfront and a 
cutoff' position r = Jfront + ^/o as a power law from the density 
in the diff'use region close to the star, up to a maximum density 
^plateau g^y^j^^ cutoff' position, the density is assumed to be 
constant, equal to n^^^^^^\ 



plateau 



do 



plateau 



if dfront < r < do 
if r > do, 

plateaux 



(A.5) 
(A.6) 



where the gas density in the plateau (n^ ) was determined by 
molecular observations (cf. references in Table A.l) and Jfront 
was evaluated at half of the peak intensity of the mid-IR emis- 
sion (Figs. 4 and 9). G^^^^ has been evaluated as explained in 
Sect. 2. The free parameters of the fit are n^, a and do, which 
define the gas density profile, and the cloud physical length L. 
The fit results for each one of the PDRs with spatially resolved 
mid-IR populations are summarised in Table A. 1 , and the mod- 
elled profiles are reported in the middle panels of Fig. 10. Along 
with n^, a and do, the model enables one to calculate the values 
of the UV radiation field Go along the cuts (Fig. 10). Good fits 
to lo are obtained (Fig. 10), except for the deepest parts of the 
PDRs, mainly at J > Jo where the density is poorly constrained. 
Since we concentrated our analysis around the AIB peak, this 
discrepancy does not aff'ect the rest of our analysis. 
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Fig. 9. Same as figure 4 but for the other sources considered in this work: from top to bottom, NGC 7023 S, NGC 7023 E, 
NGC 2023 N, NGC 2023 S, and the p-Ophiucus Eilament. The solid lines represent the cuts studied in this work, whose analy- 
sis is shown in Eig. 10. 
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Fig. 10. Compilation of our results along each PDR cut defined in Figs. 4 and 9. Upper panels: the local density and UV radiation 
field profiles as a function of the distance from the star along the r or Ay axis of Fig. A.l. Central panels: Blue, yellow, and red 
represent the contribution of PAH^, PAH^ and eVSGs respectively, as extracted from the fit. Solid black lines represent the corrected 
mid-IR emission profile. The dotted lines represent the fit obtained with the geometrical model described in Appendix A. Lower 
panels: the variation in the fraction of carbon atoms locked in eVSGs relative to the total carbon in all the AIB carriers. The 
grey-shaded regions are those that were not considered in Fig. 6. and Eq. 5. 
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